Casimir micro-sphere diclusters and three-body effects in fluids 
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Our previous article [Phys. Rev. Lett. 104, 060401 (2010)] predicted that Casimir forces induced 
by the material-dispersion properties of certain dielectrics can give rise to stable configurations of 
objects. This phenomenon was illustrated via a dicluster configuration of no n- touching objects 
consisting of two spheres immersed in a fluid and suspended against gravity above a plate. Here, we 
examine these predictions from the perspective of a practical experiment and consider the influence 
of non-additive, three-body, and nonzero-temperature effects on the stability of the two spheres. 
We conclude that the presence of Brownian motion reduces the set of experimentally realizable 
silicon/teflon spherical diclusters to those consisting of layered micro-spheres, such as the hollow- 
core (spherical shells) considered here. 



PACS numbers: 

I. INTRODUCTION 

In this paper, we investigate the influence of non- 
additive/three-body and nonzero-temperature effects on 
our earlier prediction that the Casimir force (which arises 
from quantum electrodynamic fluctuations [IH3]) can en- 
able dielectric objects (micro-spheres) with certain ma- 
terial dispersions to form stable non-touching configura- 
tions (diclusters) in fluids [4j [5] . Such micro-sphere in- 
teractions are predicted to possess a variety of unusual 
Casimir effects, including repulsive forces, [6-8 a strong 
interplay with material dispersion [4 , and strong tem- 
perature dependences [9 , and may have applications in 
microfluidic particle suspensions [lOj [11] . A typical sit- 
uation considered in this paper is depicted in Fig. [TJ 
consisting of silicon and teflon micro-spheres suspended 
in ethanol above a gold substrate. Although our ear- 
lier work considered pairs of micro-spheres suspended 
above a substrate in the additive/pair wise approxima- 
tion, summing the exact two-body sphere-sphere and 
sphere-substrate interactions, in this paper we perform 
exact three-body calculations. In Sec. [TTJ we explic- 
itly demonstrate the breakdown of the pairwise approx- 
imation for sufficiently small spheres, in which an ad- 
jacent substrate modifies the equilibrium sphere separa- 
tion, but we also identify experimentally relevant regimes 
in which pairwise approximations (and even a parallel- 
plate proximity- force/PFA approximation [12]) are valid. 
In Sec. |III[ we also consider temperature corrections to 
the Casimir interactions. Although a careful choice of 
materials can lead to a large temperature dependence 
stemming from the thermal change in the photon dis- 
tribution [9j [I3j [14], we find that such thermal-photon 
effects are negligible (< 2%) for the materials considered 
here. However, we show that substantial modifications 
to the objects separations occur due to Brownian motion 
of the micro-spheres. This effect can be reduced by low- 
ering the temperature, limited by the freezing point of 
ethanol (T « 159 K ), or by increasing the sphere diam- 
eters. We propose experimentally accessible geometries 
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FIG. 1: Schematic of two-sphere dicluster geometry consisting 
of two dielectric spheres of radii Ri and R2 separated by a 
center-center distance d from each other, and suspended by 
heights hi and /12, respectively, above a dielectric plate. 



consisting of hollow micro-spheres (which can be fabri- 
cated by standard methods [15]) whose dimensions are 
chosen to exhibit a clear stable non-touching equilibrium 
in the presence of Brownian fluctuations. We believe that 
this work is a stepping stone to direct experimental ob- 
servation of these effects. 

In fluid-separated geometries the Casimir force can be 
repulsive, leading to experimental wetting effects [TSl - 
[T8] and even recent direct measurements of the repul- 
sive force in fluids for sphere-plate geometries [T9H2T] . In 
particular, for two dielectric/metallic materials with per- 
mittivity S\ and S3 separated by a fluid with permittivity 
£2, the Casimir force is repulsive when e\ < £2 < £3 ©• 
More precisely, the permittivities depend on frequency 
lj, and the sign of the force is determined by the order- 
ing of the £fc(ift) values at imaginary frequencies uj = in 
(where Sk is purely real and positive for any causal pas- 
sive material |6]). If the ordering changes for different 
values of ft, then there are competing repulsive and at- 
tractive contributions to the force. At larger or smaller 
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separations, smaller or larger values of ft, respectively, 
dominate the contributions to the total force, and so 
the force can change sign with separation. For exam- 
ple, if s\ < 82 < £3 for large ft and e\ < £3 < £2 for 
small ft, then the force may be repulsive for small sep- 
arations and attractive for large separations, leading to 
a stable equilibrium at an intermediate nonzero separa- 
tion. Alternatively, for a sphere-plate geometry in which 
the sphere is pulled downwards by gravity, a purely re- 
pulsive Casimir force (which dominates at small separa- 
tions) will also lead to a stable suspension. These ba- 
sic ideas were exploited in our previous work [4 to de- 
sign sphere-sphere and sphere-plate geometries exhibit- 
ing a stable non-touching configuration. The effects of 
material dispersion are further modified by an interplay 
with geometric effects (which set additional length-scales 
beyond that of the separation), as well as by nonzero- 
temperature effects which set a Matsubara length-scale 
2irkT/h [14] that can further interact with dispersion in 
order to yield strong temperature corrections [9 . Ex- 
perimentally, stable suspensions are potentially appeal- 
ing in that one would be measuring static displacements 
rather than force between micro-scale objects. The sta- 
ble configurations may be further modified, however, by 
three-body effects in sphere-sphere-plate geometries and 
by Brownian motion of the particles within the potential 
well created by the Casimir interaction, and these effects 
are studied in detail by the present paper. 

Until the last few years, theoretical predictions of 
Casimir forces were limited to a small set of simple ge- 
ometries (mainly planar geometries) amenable to analyti- 
cal solution, but a number of computational schemes have 
recently been demonstrated that are capable of handling 
complicated (and, in principle, arbitrary) geometries and 
materials [22lf24j . Here, since the geometries considered 
in this paper consist entirely of spheres and planes, we are 
able to adapt an existing technique [24] based on Fourier- 
like ("spectral") expansions that semi- analytically ex- 
ploits the symmetries of this problem. This technique, 
formulated in terms of the scattering matrices of the ob- 
jects in a basis of spherical or plane waves, was developed 
in various forms by multiple authors [24-26], and we em- 
ploy the generalization of [24]. Although this process 
is described in detail elsewhere [24] and is reviewed for 
the specific geometries of this paper in the appendix, the 
basic idea of the calculation is as follows. The Casimir 
energy can be expressed via path integrals as an integral 
J °° logdet A(ft)dft over imaginary frequencies ft, where 
A is a "T-matrix" related to the scattering matrix of the 
system. In particular, one needs to compute the scatter- 
ing matrices relating outgoing spherical waves from each 
sphere (or planewaves from each plate) being reflected 
into outgoing spherical waves (or planewaves) from ev- 
ery other sphere (or plate), which can be expressed semi- 
analytically (as infinite series) by "translation matrices" 
that re-express a spherical wave (or planewave) with one 
origin in terms of spherical waves (or planewaves) around 
the origin of the new object [24]. This formalism is exact 



(no uncontrolled approximations) in the limit in which 
an infinite number of spherical/plane waves is consid- 
ered. To obtain a finite matrix A, the number of spheri- 
cal waves (or spherical harmonics Yi m ) is truncated to a 
finite order £. Because this expansion converges exponen- 
tially fast for spheres [24j[27], we find that i < 12 suffices 
for < 1% errors with the geometries in this paper. (Con- 
version from planewaves to spherical waves is performed 
by a semi-analytical formula [24] that involves integrals 
over all wavevectors, which was performed by a standard 
quadrature technique for semi-infinite integrals [28].) Al- 
though it is possible to differentiate logdetA analyti- 
cally to obtain a trace expression for the force [23], in 
this paper we use the simple expedient of computing the 
energy and differentiating numerically via spline inter- 
polation. Previously, Ref. [29] employed the same for- 
malism in order to study a related geometry consisting 
of vacuum-separated perfect-metal spheres adjacent to a 
perfect-metal plate, where it was possible to employ the 
method of images to reduce the computational complex- 
ity dramatically. That work found a three-body phe- 
nomenon in which the presence of a metallic plate re- 
sulted on a stronger attractive interaction between the 
spheres, and that this effect becomes more prominent at 
larger separations [29], related to an earlier three-body 
effect predicted for cylindrical shapes [30, 31 . Here, we 
examine dielectric spheres and plate immersed in a fluid 
and therefore cannot exploit the method of images for 
simplifying the calculation, which makes the calculation 
much more expensive because of the many oscillatory 
integrals that must be performed in order to convert be- 
tween planewaves (scattering off of the plate) and spher- 
ical waves (see appendix). We also obtain three-body ef- 
fects, in this case on the equilibrium separation distance, 
but find that the magnitude and sign of these effects de- 
pends strongly on the parameters of the problem. 



II. THREE-BODY EFFECTS 

To quantify the strength of three-body effects in the 
sphere-sphere-plate system of Fig. [I] we begin by com- 
puting how the zero-temperature equilibrium sphere- 
sphere separation d varies as a function of the sphere- 
plate separation h for two equal-radius spheres, as plotted 
in Fig. [2] To start with, we consider very small spheres, 
with radius R = 25 nm, for which the three-body effects 
are substantial. The separation dh at a given h is normal- 
ized by d^ (d as h — >• 00, i.e. in the absence of the plate). 
Several different material combinations are shown (where 
X-Y-Z denotes spheres of materials X and Y and a plate 
of material Z): polystyrene (PS), teflon (Tef), and sili- 
con (Si) spheres with gold (Au), teflon (Tef), and vacuum 
(air) plates (the latter corresponding to a fluid-gas inter- 
face). Depending on the material combinations, we find 
that dh can either increase or decrease by as much as 15% 
as the plate is brought into proximity with the spheres 
from h = 00 to h w R. (We expect even larger devia- 
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FIG. 2: Equilibrium separation d e (h)/d e (oo) between two 
R = 25nm spheres suspended in ethanol as a function of their 
surface-surface separation h from a plate (and normalized by 
the equilibrium separation for the case of two isolated spheres, 
i.e. h = oo). d e is plotted for various material combinations, 
denoted by the designation sphere-sphere-plate, e.g. a PS 
and silicon sphere suspended above a gold plate is denoted as 
PS-Si-Au. Solid/dashed lines correspond to stable/unstable 
equilibria. (In the case of a gold plate, the spheres are chosen 
to have R — 50nm.) The inset shows d e (in units of nm) for 
the case of two PS and silicon spheres (R = 57nm) above a 
gold plate. 



tions when h < R, but small separations are challenging 
for this computational method [24] and our results for 
h > R suffice here to characterize the general influence 
of three-body effects.) 

Interestingly, depending on the material combination, 
the dh can either increase or decrease as a function of 
h: that is, the proximity of the plate can either in- 
crease or decrease the effective repulsion. This is quali- 
tatively similar to previous results for vacuum-separated 
perfect-metal spheres/plates [29] in the following sense. 
Previously, the attractive interaction between a sphere 
and a plate was in general found to enhance the at- 
traction between two identical spheres as the plate be- 
came closer [29]. (There are certain regimes, not present 
here, where the attractive interaction decreases) Here, 
we observe that the sphere-plate interaction changes the 
sphere-sphere interaction with the same sign as h be- 
comes smaller: if the sphere-plate interaction is repul- 
sive, the sphere-sphere interaction becomes more repul- 
sive (larger d), and vice- versa for an attractive sphere- 
plate interaction. Since the spheres are not identical, the 
three-body effect is dominated by the sign of the stronger 
sphere-plate interaction out of the two spheres. Thus, 
examining the signs and magnitudes of the pairwise in- 
teractions in all cases of Fig.[2]turns out to be sufficient to 
predict the sign of the three-body interaction, although 
we have no proof that this is a general rule. (In contrast, 
for non-spherical objects such as cylinders, there can be 
competing three-body effects that make the sign more 
difficult to predict, even in vacuum-separated geometries 
where all pairwise interactions are attractive, which can 
even lead to a non-monotonic effect [30] I3T].) 



FIG. 3: Equilibrium separation d e (oo) (units of nm) between 
a Si sphere and either a teflon(Tef) or polystyrene(PS) sphere 
immersed in ethanol as a function of their equivalent radii R. 
Solid/dashed lines denote stable/unstable equilibria. 



Figure [2] also exhibits the interesting phenomenon of 
bifurcations, in which stable equilibria (solid lines) and 
unstable equilibria (dashed lines) appear/disappear at 
some critical h for certain materials and geometries, 
which is discussed in more detail in Sec. Ill Al As the 
sphere radius R increases, all of these three-body effects 
rapidly decrease, eventually entering an additive regime 
in which three-body effects are negligible and in which 
a parallel-plate/PFA approximation eventually becomes 



valid, as described in Sec. II B 



A. Bifurcations 

In the case of PS and Si spheres suspended above ei- 
ther a gold or teflon plate, one can observe the emer- 
gence or disappearance of a stable (solid) and unstable 
(dashed) pair of equilibria as h decreases from h = oo, 
respectively, as evidenced by the blue curves in Fig. [2] 
(teflon plate) and Fig. [2] (inset) (Au plate). This can 
be qualitatively explained by the fact that the isolated 
sphere-sphere interactions exhibit a natural bifurcation 
for sufficiently- large spheres, in conjunction with the fact 
that the presence of the plate typically acts to either in- 
crease or decrease the sphere-sphere interaction, depend- 
ing on the sign of the dominant sphere-plate interaction, 
as explained above. 

In particular, Fig. [3] shows the isolated Si-PS and Tef- 
Si sphere-sphere equilibrium separation d e as a function 
of the radius R of the spheres. As a consequence of 
its material dispersion (similar to phenomena observed 
in U), the Si-PS combination exhibits a bifurcation at 
R « 55 nm where the stable and unstable equilibria, 
such that there is no equilibrium for larger R (the in- 
teraction is purely attractive). The Tef-Si combination 
exhibits no such bifurcation (even if we extend the plot 
to R = 300 nm), because it has no unstable equilibrium: 
the interaction is purely repulsive for small separations 
and attractive for large separations. Therefore, if the Si- 
PS radius is above or below the 55 nm bifurcation, the 
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presence of the plate can shift this bifurcation and lead 
to a bifurcation as a function of h as in Fig. |2j whereas 
no such bifurcation with h appears for Tef-Si. 

In the Si-PS-Au case of a gold plate with Si-PS 
spheres, the sphere-plate interactions turns out to be 
primarily repulsive, which should push the bifurcation 
in Fig. [3] to the right (shrinking the attractive region) 
as h decreases. Correspondingly, if we choose a radius 
R = 57 nm just to the right of isolated-sphere bifur- 
cation, then as h decreases the Si-PS-Au combination 
should push the bifurcation past R = 57 nm leading to 
the creation of a stable/unstable pair for small h, and 
precisely this behavior is observed in the inset of Fig. |2j 
Conversely, for the Si-PS-Tef case of a teflon plate with 
Si-PS spheres, the sphere-plate interaction is primarily 
attractive, and the opposite behavior occurs: choosing a 
radius R = 25 nm to the left of the isolated-sphere bifur- 
cation, decreasing h increases the attraction and moves 
the bifurcation to the left in Fig. [3] eventually causing 
the disappearance of the stable/unstable equilibrium at 
R = 25 nm. Correspondingly, for the Si-PS-Tef curve 
in Fig. [2j we see the disappearance of a stable/unstable 
pair for sufficiently small h. 
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FIG. 4: Center-surface L e (solid lines) and surface-surface h e 
(dotted lines) equilibrium separation (units of nm) of a teflon 
(red) or Si hollowed sphere (shown on the inset) suspended in 
ethanol above a gold plate, as a function of radius R (units of 
nm) . The equilibria are plotted for different values of the fill- 
fraction a, denned as the ratio of the spherical-shell thickness 
over the radius of the sphere. Solid/dashed lines correspond 
to stable/unstable equilibria. 



B. The additive regime 

In general, three-body effects can expected to disap- 
pear in various regimes where key parameters of the in- 
teraction become small. First, for large radii, where h 
(the sphere-plate separation) and d (the sphere-sphere 
separation) become small compared to i?, eventually the 
Casimir interaction is dominated by nearest-surface in- 
teractions, or the proximity-force approximation (PFA), 
in which the force can be approximated by additive 
surface-surface "parallel-plate" forces [12, 32 , 33 . In or- 
der to damp the Brownian fluctuations as described in 
the next section, we actually propose to use much larger 
{R > 5 jam) spheres, and we quantify the accuracy of 
PFA in this regime below. Second, as h becomes large 
compared to d, the effect of the plate becomes negligi- 
ble and three-body effects disappear; this is apparent in 
Fig. [2j where d e — >• d e (oo) when h ^> d e . Third, in the 
limit where one of the spheres is much smaller than the 
other sphere, then the smaller sphere has a negligible ef- 
fect on the sphere-plate interaction of the larger sphere, 
and at least some of the three-body effect disappear as 
described below. In fact, we find that even for a situa- 
tion in which one sphere is only a few times smaller than 
the other, the three-body effects tend to be negligible. 
For the sphere-radius regime considered in our previous 
work, we argue below that equal-height suspension of the 
two spheres leads to a strong asymmetry in sphere radii 
that tends to eliminate three-body effects. 

To begin with, let us consider sphere radii on the order 
of 10 2 nm, as in our previous work [4] . We wish to make a 
bound dicluster, at some separation d, of two spheres (Si 
and teflon) that are suspended above a gold substrate 



by Casimir repulsion in balance with gravity. Further- 
more, suppose that we wish to suspend both spheres at 
the same equilibrium height /i e , and therefore choose the 
radii of the two spheres to equate their h e values. In 
Fig. |1J we plot h e as a function of radius R for the iso- 
lated sphere-plate geometries (d — >• 00). For example, 
with an Si sphere of radius R = 100 nm, the (stable) 
equilibrium height is h e = 298.17 nm, whereas to ob- 
tain the same h e value for teflon one needs a much larger 
teflon sphere of radius R = 217.2 nm, primarily because 
the Casimir repulsion is stronger for teflon. If, instead of 
a pairwise calculation, we perform an exact three-body 
calculation of the h e values for these radii at the equi- 
librium sphere-sphere separation d e = 92.8 nm, we find 
that the h e values change by < 1%. Conversely, if we 
keep h e fixed and compute the three-body change in d e 
(compared to h — >• 00), again we find that the change 
is < 1%. As mentioned above, the small size of the 
Si sphere makes it unsurprising that the Si sphere does 
not change the equilibrium h e of the much larger teflon 
sphere. Furthermore, the sensitivity of the sphere-sphere 
force Fd to the teflon h is equal to the sensitivity of the 
teflon sphere-plate force to d, thanks to the equiva- 
lence dF d /dh = -d 2 U/dddh = dF h /dd where U is the 
energy. Therefore, one would also not expect the finite 
value of h e for the Si sphere to modify the equilibrium 
d e . Size asymmetry alone, however, does not explain 
why the finite h e of the teflon sphere does not affect the 
sphere-plate interactions of the Si sphere. Even if the Si 
sphere were of infinitesimal radius, the Casimir-Polder 
energy the Si sphere would be determined by a Green's 
function at the Si location [2], and if the Si sphere is 
at comparable distance d e ~ h e from both the teflon 
plate and the sphere, one would in general expect the 



response to a point-dipole source at the Si location (the 
Green's function) to depend non-additively on the teflon 
sphere and the plate even for an infinitesimal Si sphere. 
However, in the present case we do not observe any non- 
additive effect on the Si-sphere h e , because the factor 
of three (approximately) difference between h e and d e is 
already sufficient to eliminate three-body effects (as in 
Fig.§. 

Figure [4] also exhibits a bifurcation of stable (solid 
lines) and unstable (dashed lines) equilibria that causes 
the stable h equilibrium to vanish for Si at large radii. 
In order to utilize larger spheres to reduce the effects of 
Brownian motion in the next section, one can consider in- 
stead a geometry of hollow air-filled spherical shells with 
outer radius R and shell thickness aR (so that a = 1 
gives a solid sphere) . Such hollow microspheres are read- 
ily fabricated with a variety of materials [15]. As Fig. [4] 
shows, decreasing the shell thickness a pushes the bi- 
furcation to larger R, and also increase h e by making 
the sphere more buoyant. This modification allows us 
to consider R « 10 /am in the next section, where PFA 
should be accurate. For only 3 fim spheres and 500 nm 
separations in fluids, we previously found that the cor- 
rection to PFA (which scales as d/R to lowest order [34l - 
[36] ) was only about 15%. For the three times larger radii 
and somewhat smaller separations in the next section the 
corrections to PFA are typically < 5%, sufficient for our 
current purposes. 



III. NONZERO TEMPERATURE AND 
EXPERIMENTS 



In this section, we address a number of questions 
of consequence to an experimental realization of the 
teflon/silicon two-sphere dicluster of Fig. [I] In particu- 
lar, we consider several ways in which a nonzero temper- 
ature can disrupt the observation of stable equilibria. A 
nonzero temperature will manifest itself in at least two 
important ways. First, there will be a change in the 
Casimir force between the objects due to the presence 
of real (non-virtual) photons in the system. Second, the 
inclusion of nonzero temperature will cause the spheres 
to experience Brownian motion arising from the thermal 
agitations in the fluid [37 . We consider the influence of 
both of these effects on the observability of stable particle 
clusters and suspensions. 

At zero temperature, the Casimir force F is determined 
by an integral F — J °° of a complicated integrand 

/(£) evaluated at imaginary frequencies £ [2]. At T > 0, 
the integral is replaced by a finite sum over Matsubara 
frequencies u n = 27rnkT/H, arising from the poles of the 
coth photon distribution along the imaginary frequency 
axis p~4j[38], leading to a force Ft given by: 
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FIG. 5: Geometry of a hollow air-filled core sphere suspended 
above a layered plate with layer thickness H. Explicitly shown 
is the thickness dimension as a function of a 



which is exactly a trapezoidal-rule approximation to the 
zero-temperature force with a discretization error deter- 
mined by the Matsubara wavelength At = 27tc/£t = 
H/kT [39]. Because the integrand /(£) is smooth and 
typically varies on a scale much slower than 1 / At , where 
At = 7.6 /im at room temperature T = 300 K, the finite- 
T correction to the zero-temperature Casimir force is 
often negligible [3]. However, in fluids, as is the case 
here, larger temperature effects have been obtained [9 
by dispersion-induced oscillations in /(£), and so we must 
check our previous zero-temperature predictions against 
finite-T calculations. For the Tef-Si-Substrate case con- 
sidered here, we find that T > corrections to the T = 
forces are no more than 2% over the entire range of sepa- 
rations considered here, and hence they can be neglected. 

The presence of Brownian motion proves a much more 
difficult experimental complication to overcome. First, 
Brownian motion will lead to random fluctuations in the 
position of the spheres, making it hard to measure their 
stable separations in an experiment [37 . Second, and 
more importantly, sufficiently large fluctuations can drive 
the Si sphere to "tunnel" past its unstable equilibrium 
position with the gold plate, leading to stiction [37] since 
the Si-Au interaction is purely attractive for small sepa- 
rations. The remainder of this section will revolve around 
the question of how and whether one can overcome both 
of these difficulties to observe suspension in experiments. 
In particular, we consider observation of the average sep- 
aration of the spheres over a sufficiently long time, but 
not so long that stiction occurs, and analyze the separa- 
tion statistics and the stiction timescale. First, however, 
we describe how the parameters are chosen so that Brow- 
nian fluctuations are not so severe. 

The sphere geometry that we consider is depicted in 
Fig. [5] a hollow spherical shell suspended by a surface- 
surface separation h above a layered substrate, consisting 
of a thin indium tin oxide (ITO) film of thickness H de- 
posited on a gold substrate, where the purpose of the 
ITO layer is to eliminate the Si-sphere instability as ex- 
plained below. The thickness of the shell is denoted as 
t = aR, where a is a convenient fill-fraction parame- 
ter. We consider hollow spheres in order to increase R 
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and thereby reduce Brownian fluctuations. In particu- 
lar, both the Brownian fluctuations and the probability 
of stiction in the case of the silicon sphere are reduced 
by increasing the strength of the Casimir force, which 
can be achieved by increasing R since the Casimir force 
scales roughly with surface area, and below we consider 
radii from 1 to 10 jam. In this regime, as quantified in 
the previous section, a simple PFA approximation is suf- 
ficient to accurately compute the forces and separations. 
However, because the gravitational force scales as i? 3 , for 
large R the gravitational force will overcome the Casimir 
force and push the Si sphere past its unstable equilib- 
rium into stiction. In order to reduce the gravitational 
force while keeping the surface area fixed, we propose 
using a hollow Si sphere. We find that in addition to 
hollowing the spheres, it is also beneficial to deposit a 
thin ITO film,o n top of the gold substrate (the permit- 
tivity of ITO is modeled via an empirical Drude model 
with plasma frequency uj p = 1.46 x 10 15 rad/s and decay 
rate 7 = 1.53 x 10 14 rad/s). The ITO layer acts to de- 
crease the equilibria separations and therefore increase 
the Casimir interactions between the spheres and the 
substrate. However, because the Casimir force between 
teflon/silicon and ITO is attractive/repulsive at small 
separations, respectively, increasing H pushes the silicon- 
substrate's unstable equilibrium to smaller separations 
while introducing a teflon-substrate unstable equilibrium 
that gets pushed to larger separations. In what follows, 
we find that H from 14-30 nm is sufficient to obtain ex- 
perimentally feasible suspensions, although here we only 
consider the case of H = 15 nm. 

The effect of hollowing the spheres is shown in the top 
panel of Fig. [6] smaller a values push the stable/unstable 
bifurcation of teflon out to larger R. Hollowing the sil- 
icon sphere is not necessary because silicon has no un- 
stable equilibrium (in this configuration it is repulsive 
down to zero separation). However, as shown in the bot- 
tom panel of Fig. [6j hollowing the silicon sphere does 
change its h e at a given R. For example, one can choose 
a Tef a = 0.142 and a Si (a = 0.14) to obtain the same 
equilibrium surface-to-center height L e over a wide range 
of sphere radii, as shown in upper-right inset of Fig. [6] 
Alternatively, one can choose a hollow teflon sphere to 
match the equilibrium surface-surface separations h e for 
equal sphere radii, as shown in the lower-right inset of 
Fig.© 



A. Statistics of Brownian motion 

As mentioned above, Brownian motion will disturb the 
spheres by causing them to move randomly about their 
stable equilibrium positions, and this can cause the Si 
sphere to move past its unstable equilibrium point, in- 
ducing it to stick top the plate. To quantify the range 
of motion of both spheres about their equilibria, we con- 
sider the statistical properties of their fluctuations. In 
particular, we consider the average plate-sphere separa- 
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FIG. 6: Surface-surface equilibrium height h e (units of nm) 
for the hollowed-sphere geometry of Fig. [5] consisting of either 
a Si (top) or teflon (bottom) hollowed sphere (fill-fraction 
a) suspended in ethanol above a H = 15 nm ITO layered 
gold plate, as a function of sphere radius R (in units of /mi). 
Solid/dashed lines correspond to stable/unstable equilibria. 
h e is plotted for different values of a, denoted in the figure. 
The top inset plots the center-surface separation L e (in units 
of nm) as a function of R of a hollowed teflon/Si (red/blue 
lines) sphere suspended again above a gold plate, for a = 
0.14/0.142. The lower inset shows h e for both teflon and Si 
spheres for R £ [8.6, 10]/im. 



tions (K)t and average sphere-sphere separations (oI)t 
near room temperature (T = 300 K), determined by an 
ensemble average over a Boltzman distribution. For ex- 
ample, (h)r is given by: 



J °° dzzexp(U(z)/kT) 
J °° dz exp (U(z)/ kT) ' 



(2) 



where U (z) is the total energy (gravity included) of the 
sphere-plate system at a surface-surface height z. (A 
similar expression yields (c2)t)- I n the case of teflon, the 
short-range attraction means that the suspension is only 
metastable under fluctuations; here, we only average over 
separations prior to stiction by restriction z to be > the 
unstable equilibrium, and consider the stiction timescale 
separately below. In addition to the average equilibrium 
separations, we are also interested in quantifying the ex- 
tent of the fluctuations of the spheres, which we do here 
by computing the 95% confidence interval {<j_,<j + }, de- 
fined as the spatial region over which the sphere is found 
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FIG. 7: Average (/i) (thick lines) and equilibrium h e (thin 
lines) height (in units of nm) of a hollowed teflon/Si (blue/red 
lines) sphere suspended above a H = 15 nm ITO layered 
gold plate, for a = 0.142/0.13, as a function of sphere ra- 
dius R (in units of /xm). Solid/dashed lines correspond to 
stable/unstable equilibria. The red/blue shaded regions in- 
dicate positions where the teflon/Si spheres are found with 
95% probability. The inset shows (d) (thick line) and d e (thin 
line) separations as a function of their radius for two equal 
radii teflon Si spheres. The gray shaded region indicate the 
separations which the teflon and Si spheres are found with 
95% probability. 



with 95% probability around the equilibria, where <j± de- 
notes the lower/upper bound of that interval. These re- 
sults are shown in Fig.[7]for h, with d shown in the inset, 
in which shaded regions indicate the confidence intervals, 
as a function of R where a is chosen to yield approxi- 
mately equal h e (a = 0.142 for teflon and a = 0.13 for 
Si). (Note that the horizontal separation (d) is a purely 
Casimir interaction and the difference here from a = 1 is 
negligible in the PFA regime.) As predicted above, the 
Brownian fluctuations of the spheres vanish as R —> oo 
and are dramatically suppressed for R > 5 /im, where one 
finds (h) « h e . In addition, we find that the teflon sphere 
can safely avoid the unstable equilibrium and stiction in 
the sense that the unstable equilibrium is far outside the 
confidence interval; the timescale of the stiction process 
is quantified below. The asymmetrical nature of the con- 
fidence interval results from the fact that the Casimir 
energy decreases as a function of z, and as a consequence 
the Brownian excursions favor the +z direction. The 
fluctuations in (d) are substantially larger than those in 
(h) (nor is there any obvious reason why they should be 
comparable, given that the nature of the sphere-sphere 
equilibrium is completely different from the sphere-plate 
equilibrium), making the precise value of d e potentially 
harder to observe. 

Instead of considering the Brownian statistics as a 
function of i?, we can instead consider the statistics as 
a function of a for fixed radii « 10/im (chosen to ob- 
tain nearly equal sphere-center heights L e ), as shown in 
Fig. [8] One key point is that there is a minimum allowed 
a: if a is too small, the buoyant force (assuming an air- 
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FIG. 8: Average (h) (thick line) and equilibrium h e (thin 
line) height (in units of nm) of a hollowed Si/teflon sphere 
(top/bottom) of radii R = 10/9.915/xm suspended in ethanol 
above a H = 15 nm ITO layered gold plate, as a function of 
fill- fraction a (indicated in Fig. [SJ. Shaded regions indicate 
h positions where the Si/teflon spheres are found with 95% 
probability. Solid/dashed lines indicate stable/unstable equi- 
libria. For reference, we state the equilibrium d e and average 
(d) horizontal seperations between R — 10/9.915/xm Si/Tef 
spheres in the top figure. 



filled hollow sphere) will eventually become positive and 
the sphere will float, although this limitation is removed 
if one could infiltrate the hollow sphere with the fluid. 
For the teflon sphere, there is also an upper limit to a 
for a given R to avoid stiction as discussed previously. 



B. Stiction and tunnelling rates 

As mentioned above, the stable equilibrium for the 
teflon sphere is actually only metastable — because the 
Casimir force is attractive for small separations, given a 
sufficiently long observation time r the sphere will "tun- 
nel" (via Brownian fluctuations) past the energy barrier 
A posed by the unstable equilibrium, and stick to the 
plate (stiction). Given the energy barrier, the tempera- 
ture T, and the viscous drag on the particle, we can apply 
standard methods [37, 40, 41 to compute the timescale 
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FIG. 9: Energy barrier A/kT of a hollowed teflon sphere 
suspended in ethanol above a H — 15 nm ITO layered gold 
plate at T — 300 K, as a function of sphere radius R (in 
units of /im) and for different values of fill- fraction a. The 
inset shows the energy landscape U/kT as a function of the 
surface-surface height h (units of nm) for a teflon sphere of 
radius R = 10/xra with a fill fraction of a = 0.142. 



for stiction. This calculation, which is described in detail 
below, shows that for various values of the fill factor a 
the expected time r to stiction (which increases exponen- 
tially with A/kT) can vary dramatically, but can easily 
be made on the order of years. 

The energy barrier A/kT is plotted versus the teflon 
sphere radius R for various a in Fig. [9j and can easily 
be made > 10 to obtain a very long metastable lifetime. 
As we discussed earlier, the A increases with R at first 
because this increases the Casimir force, but has a max- 
imum at some R where gravity begins to dominate. De- 
creasing a decreases the gravitational force and therefore 
increases both the maximum A and the corresponding R. 
A typical energy landscape U(z) is shown in the inset, ex- 
hibiting a local minimum at a height h e and an unstable 
equilibrium (maximum) at h u . Also noted on the inset is 
the "tunneling" height h* > h e at which U(h*) = U(h u ). 
Figure [9] also shows the energy barrier A/kT of a sil- 
icon sphere (a = 0.1288 « a c , R = 10 /xm) in the ab- 
sence of the ITO layer (H = 0) to be significantly smaller 
than that of teflon. Of course A/kT in this case could 
be made larger merely by choosing a « a c , but we find 
(below) that achieving experimentally realizable lifetimes 
severely limits the range of realizable a, i.e. requires that 
the Si thickness be known to within a few nanometers. 

Because A ^> fcT, the lifetime r of a Brownian particle 
trapped around a local minimum of a potential U (z) can 
be approximated by [40] : 



a A/kT 



7 \l/2 

4^) 



1_ 
2uj 



27T fjS 



(3) 



where 7 is the viscous drag coefficient (drag force = 
—7 velocity), lj and ft characterize the curvature of U(z) 
at the energy maximum and minimum respectively [as 
defined in Eq. d5b], C(6) is a transcendental function de- 
fined in Eq. (pJ7 and S is an integral of the potential 



barrier defined by Eq. Q. Let m be the mass of the 
sphere. The drag coefficient for a sphere of radius R in a 
fluid with viscosity 77 is 7 = 67rRr]/m [42], where a typical 
viscosity is r] « 1.17 ± 0.06 mPas for ethanol [43 . The 
other quantities are given by: 



5 = 2 / dzy/-2mU(z) 



U"(h u ) 



m 



U"{h e ) 
m 



£(5) = exp 



2 ^l 2 

-J dzln(l-e- 6/Acos2z ) 



(4) 

(5) 
(6) 



Combining these formulas and choosing different values 
of R and a to obtain different barriers A and landscapes 
U (z) as in Fig. [9J the lifetime r can be designed to take 
on a wide range of values. The exponential dependence 
on A means that r rapidly transitions from very short 
to very long as a changes, but can easily be made large. 
For example, with R = 8.5/im and a < 0.15, one ob- 
tains r > 40 days. (Conversely, for sufficiently large a 
one could design experiments where stiction occurs on an 
arbitrarily fast timescale, but in this A ~ kT regime the 
approximations of Eq. (|3| are no longer valid.) 

Strictly speaking, this is a conservative estimate of the 
timescale because the drag coefficient 7 for a sphere above 
a plate is larger than that of an isolated sphere. As the 
sphere approaches the plate, the drag is dominated by 
the "lubrication" problem of the fluid squeezed between 
the sphere and the plate, and the drag increases dramat- 
ically [44]. 



IV. CONCLUSION 

Even including the thermal motion of the particles and 
the finite lifetime of metastable suspensions, the stable 
suspension and separation of particle diclusters appears 
to be experimentally feasible. In the experimentally rel- 
evant regimes, these effects consist primarily of pairwise 
sphere-sphere and sphere-plate interactions; while three- 
body effects become significant for smaller spheres, the 
increased Brownian fluctuations for small spheres makes 
such an experiment challenging. Although the systems 
considered here consisted of silicon and teflon spheres 
above layered substrate in ethanol, many other materi- 
als combinations could potentially be explored to modify 
these phenomena, including multi-material sphere sys- 
tems such as multi-layer spheres or patterned substrates 
that could exhibit unusual effective dispersion phenom- 
ena. Although we considered hollow (air core) spheres, 
one could also use fluid-filled spheres or similar modifi- 
cations in order to modify the effect of gravity. Alter- 
natively, one could use non-spherical geometries such as 
disks, which have a both surface area and volume propor- 
tional to R 2 so that gravity does not dominate asymp- 
totically. We have recently demonstrated computational 
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methods capable of accurate modeling of such geome- 
tries, and find that the additional rotational degrees of 
freedom can lead to additional phenomena such as tran- 
sitions in the stable orientation with separation [45]. In 
general, the possibility of both repulsion and stable equi- 
libria in fluids (whereas the latter are not possible in vac- 
uum [46 but do exist in the critical casimir fluids [471 148]) 
opens the possibility of a rich and currently little explored 
territory for Casimir physics, and it is likely that many 
effects remain to be discovered. 



Appendix 

In what follows, we write down an expression for the 
Casimir energy of of the system in Fig. [I] in terms of 
the scattering and translation matrices of the individual 
objects (spheres and plates) of the geometry. A simi- 
lar expression was derived in [29] in the case of perfect- 
metal vacuum-separated objects, for which an additional 
simplification, based on the method of images, was pos- 
sible [49 . Here, we consider the more general case of 
fluid-separated dielectric objects. 

The starting point of the Casimir-energy expression 
is the well-known scattering-matrix formalism, derived 
in [24] [25], in which the Casimir energy U between an 
arbitrary set of objects can be written as: 



U 



__ h ° r 



dn log det 



-l 

OO 5 



(7) 



where M^ 1 = diag(Fi,F 2 , ...) and the matrix M is given 
by: 



Ff 1 X 12 X 13 
X 21 F^ 1 X 23 



(8) 



where F^(ft) is the matrix of inside/outside scattering am- 
plitudes of the zth object, and the translation matrix 
that relates the scattering matrix of the ith and jth ob- 
jects, as described in [24 J. Here, the plate is labeled by 
the index i = 1 whereas the left and right spheres are 
labeled as i — 2 and i = 3, respectively. 

For computational convenience, the determinant in 
Eq. can be re-expressed in terms of standard oper- 
ations on the block matrices composing M, and in this 
case we find that: 



det 



where 



f oo = det (x - A/" (1) ) det (x - A/" (2) ) 

det (x - (i - x^y 1 a(x- M^y 1 B 



A/" (2) = F 3 X 31 F!X 13 , A = F 3 X 32 - F 3 X 31 F!X 12 ; 

B = F 2 X 23 - F 2 X 21 F 1 X 13 , A/" (1) = F 2 X 21 F 1 X 12 , (9) 



where (X - A/" (1) ) and (X - A/" (2) ) yield the individual 
interaction energies of the left and right spheres with the 
plate, respectively. Because of the logarithm in Eq. (|7|), 
it is possible to re-express the energy as: 



U = -^S 2 (h 2 ) + £int(hi,h 2 ,d), 



(10) 



where, 



he 
2^ 



dn log det (x-N {l) ^J 
he r°° / \ 

£ 2 (h 2 ) = ^-J etelogdet (l - M {2) ) , 



(11) 



are the individual interaction energies of the left (1) and 
right (2) spheres above a plate, in the absence of the 
other sphere, and £i n t(hi^h 2 ^d) is a three-body interac- 
tion term given by: 



he 

2^ 



dn log det 



(12) 



Finally, for completeness, we write down simplified ex- 
pressions for the intermediate matrices N^ l \ A and £>, 
in terms of appropriate and rapidly-converging multipole 
and Fourier basis, as explained in [24j. The expression 
for £i 5 2 was derived in [24] and thus here we can simply 
quote the result for the matrices and JV^ 2 \ In par- 
ticular, [24] expresses the matrices in terms of a spherical 
multipole basis, indexed by the quantum numbers Z, m, 
and P, corresponding to angular momentum, azithmutal 
angular momentum, and polarization [TE (P = E) or 
TM (P = M)\. The matrices A/" (i) are given by: 



Af {j) -S ,T K 



I 



ee(j') 
lmP,h 

00 k ± dk ± e - 2 W k i+" 2 



2tt 2*V*a + 
^2Di m p ik±Q r Q Dl ±Q l/rnfp/ (2£ Q?P / - 1) , 

(13) 



Q 



where is the Fourier momentum parallel to the plate, 



the T\ 



ImPdmP 



are the outside scattering amplitudes of 
sphere j, are the planar reflection coefficients (Fresnel 
reflection coefficients in the case of an isotropic plate), 
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and Di m p^ ±m are conversion matrices: 



Dl m E,k ± E = Di m M,k ± M 

|k±' 



1 4n(2l + 1){1 - m)\ 
l(l + l)(l + m)\ 



x 



lmM,k±E 



-D, 



lmE,k ± M = —im\ 



/47r(2Z + l)(Z-m)! 



Z(Z + l)(Z + ra)! 
x A e -*™<^ P ™ (^\s? ± + K 2 /i?j , (14) 

given in terms of associated Legendre polynomials P™ 
and their derivatives with respect to their corresponding 
argument P^ m . 

Upon a number of algebraic manipulations, similar ex- 
pressions can be obtained for the matrices A and S, not 
found in previous works, and in particular we find that: 

A — T ee 1 / 23 

SilmPJ'm'P' — J R,lmP,lmP u lmP,V m' P' 

/ ~i\m —m -m' —m -pee o 
T \ — 1) I J- RJmPJmPPlmPJ'm'P' 

(15) 
(16) 



+ i 



-ree i /32 

J~ T, lrn P,lmP U h 

^L^lmPdmpPlmP.l'm 1 P' , 



-BlmPJ'm'P' — J~ LlrnPlrnP^lrnPl'm' P 
•m' — m -ree 



where 



f°° k ± dk ± 

JO W 



-(h 2 +h 3 )^k 2 ± +K 2 



2Ky/k 2 ± + K 2 

x DimP,k ±Q r Q Di ±Q Vrnfpi (25 Q:P , - 1) , 

(17) 



Q 



and where the J m (Sk±) is a Bessel function of the first 
kind evaluated at different values of Sk± , where S is given 
by the projection of the sphere center-center separation 
onto the plate axis: 

S = ^(d + R x + R 2 ) 2 - (hi + R x - h 2 - R 2 ) 2 (18) 



From a numerical perspective, all that remains in or- 
der to obtain the Casimir energy in Eq. ^ is to evalu- 
ate the various matrix entries and perform standard nu- 
merical operations, such as inversion and multiplication, 
which we perform using standard free software [50} . For 
the small matrices that we consider, most of the time is 
spent evaluating the various matrix elements, which can 
be numerically expensive due to the integration of the 
oscillatory Bessel functions in A and S, although spe- 
cialized methods for oscillatory and Bessel integrals are 
available that may accelerate the calculation [5lJ [52] . 



[1] H. B. G. Casimir, Proc. K. Ned. Akad. Wet. 51, 793 
(1948). 

[2] E. M. Lifshitz and L. P. Pitaevskii, Statistical Physics: 

Part 2 (Pergamon, Oxford, 1980). 
[3] K. A. Milton, Journal of Physics A: Mathematical and 

General 37, R209 (2004). 
[4] A. W. Rodriguez, A. P. McCauley, D. Woolf, F. Capasso, 

J. D. Joannopoulos, and S. G. Johnson, Phys. Rev. Lett. 

104, 160402 (2010). 
[5] A. W. Rodriguez, J. Munday, D. Davlit, F. Capasso, J. D. 

Joannopoulos, and S. G. Johnson, Phys. Rev. Lett. 101, 

190404 (2008). 

[6] I. E. Dzyaloshinskii, E. M. Lifshitz, and L. P. Pitaevskii, 

Adv. Phys. 10, 165 (1961). 
[7] J. Munday, F. Capasso, and V. A. Parsegia, Nature 457, 

170 (2009). 

[8] O. Kenneth, I. Klich, A. Mann, and M. Revzen, Phys. 
Rev. Lett. 89, 033001 (2002). 



[9] A. W. Rodriguez, D. Woolf, A. P. McCauley, F. Capasso, 
J. D. Joannopoulos, and S. G. Johnson, Phys. Rev. Lett. 
105, 060401 (2010). 
[10] S. J. Rahi and S. Zaheer, arXiv:cond-mat/0909.4510vl 
(2009). 

[11] A. P. McCauley, A. W. Rodriguez, J. D. Joannopoulos, 

and S. G. Johnson, Phys. Rev. A 81, 012119 (2010). 
[12] B. V. Derjaguin, I. I. Abrikosova, and E. M. Lifshitz, Q. 

Rev. Chem. Soc. 10, 295 (1 956), URL |http: //dx.doi . | 

org/10. 1039/Q R956100 02951 
[13] M. Bostrom and B. E. Sernelius, Phys. Rev. Lett. 84, 

4757 (2000). 

[14] M. Bordag, B. Geyer, G. L. Klimchitskaya, and V. M. 
Mostepanenko, Phys. Rev. Lett. 85, 503 (2000). 

[15] D. Wilcox, M. Berg, T. Bernat, D. Kellerman, and 
J. Cochran, Hollow and Solid Spheres and Micro- 
spheres -.Science and Technology Associated with Their 
Fabrication and Application, Vol 372 (Materials Re- 



11 



search Society, 1995). 

[16] D. Bonn, J. Eggers, J. Indekeu, J. Meunier, and E. Rolley, 
Rev. Mod. Phys. 81, 739 (2009). 

[17] P.-G. d. Gennes, F. Brochard-Wyart, and D. Quere, Cap- 
illarity and Wetting Phenomena: Drops, Bubbles, Pearls, 
Waves (Springer, 2004). 

[18] J. N. Israelachvili, Intermolecular and Surface Forces (El- 
sevier, 2011). 

[19] U. Mohideen and A. Roy, Phys. Rev. Lett. 81, 4549 
(1998). 

[20] J. N. Munday and F. Capasso, Phys. Rev. A 75, 
060102(R) (2007). 

[21] A. A. Feiler, L. Bergstrom, and M. W. Rut- 
land, Langmuir 24, 2274 (2008), pMID: 18278966, 
http://pubs.acs.org/doi/pdf/10.1021/la7036907, U RL 
|http : //pubs . acs . or g/doi/abs/1 . 1021/l a7036907 

[22] A. Rodriguez, M. Ibanescu, D. Iannuzzi, J. D. 
Joannopoulos, and S. G. Johnson, Phys. Rev. A 76, 
032106 (2007). 

[23] M. T. H. Reid, A. W. Rodriguez, J. White, and S. G. 

Johnson, Phys. Rev. Lett. 103, 040401 (2009). 
[24] S. J. Rahi, T. Emig, R. L. Jaffe, and M. Kardar, Phys. 

Rev. D 80, 085021 (2009). 
[25] T. Emig, N. Graham, R. L. Jaffe, and M. Kardar, Phys. 

Rev. Lett. 99, 170403 (2007). 
[26] O. Kenneth and I. Klich, Phys. Rev. B 78, 014103 (2008). 
[27] A. Canaguier-Durand, P. A. Maia Neto, I. Cavero-Pelaez, 

A. Lambrecht, and S. Reynaud, Phys. Rev. Lett. 102, 

230404 (2009). 

[28] A. C. Genz and A. A. Malik, SIAM J. Numer. Anal. 20, 
580 (1983). 

[29] P. Rodriguez-Lopez, S.J. Rahi, and T. Emig, Phys. Rev. 

A 80, 022519 (2009). 
[30] S. J. Rahi, A. W. Rodriguez, T. Emig, R. L. Jaffe, S. G. 

Johnson, and M. Kardar, Phys. Rev. A 77, 030101 (R) 

(2008). 

[31] A. Rodriguez, M. Ibanescu, D. Iannuzzi, F. Capasso, 

J. D. Joannopoulos, and S. G. Johnson, Phys. Rev. Lett. 

99, 080401 (2007). 
[32] M. Bordag, U. Mohideen, and V. M. Mostepanenko, 

Phys. Rep. 353, 1 (2001). 
[33] H. Gies and K. Klingmiiller, Phys. Rev. Lett. 96, 220401 

(2006). 

[34] P. A. Maia Neto, A. Lambrecht, and S. Reynaud, Phys. 



Rev. A 78, 012115 (2008). 
[35] M. Francisco D., L. Fernando C., and V. Paula I, Journal 

of Physics: Conference Series 161, 012015 (2009). 
[36] H. Gies and K. Klingmuller, Phys. Rev. Lett. 96, 220401 

(2006). 

[37] H. Risken, The Fokker-Plank Equation: Methods of Solu- 
tion and Applications (Springer- Verlag, Heidelberg New 
York, 1996). 

[38] J. Schwinger, L. L. DeRaad, and K. Milton, Ann. Phys. 
115, 1 (1978). 

[39] J. P. Boyd, Chebychev and Fourier Spectral Methods 
(Dover, New York, 2001), 2nd ed. 

[40] V. I. Mel'nikov, Physics Reports 209, 1 (1991), 
ISSN 370-1573, URL [http://www. sciencedirectT 
com/science/article/B6TVP-46SXR6S-DN/2/ 
2738dl 1 ledbf 3c 995d7f 494940a32ceel 

[41] S. Chandrasekhar, Rev. Mod. Phys. 15, 1 (1943). 

[42] L. D. Landau and E. M. Lifshitz, Fluid Mechanics (Perg- 
amon Press, Elmsford, New York, 1987). 

[43] D. G. Friend and M. L. Huber, International Jour- 
nal of Thermophysics 15, 1279 (199 4), ISSN 0195- 
928X, 1 0.1007/BF014 58836, URL |http://dx. doi.org/ 
10.1007 /BF01458836 

[44] B. Hamrock, Fluid Film Lubrication (McGraw-Hill, New 
York, 1994). 

[45] M. T. H. Reid, , J. White, and S. G. Johnson, 

arXiv:1010.5539vl (2010). 
[46] S. J. Rahi, M. Kardar, and T. Emig, Phys. Rev. Lett. 

105, 070404 (2010). 
[47] M. Trondle, S. Kondrat, A. Gambassi, L. Harnau, and 

S. Dietrich, J. Chem. Phys. 133, 074702 (2010). 
[48] T. F. Mohry, A. Maciolek, and S. Dietrich, Phys. Rev. E 

81, 061117 (2010). 
[49] L. S. Brown and G. J. Maclay, Phys. Rev. 184, 1272 

(1969). 

[50] R. K. S. Hankin, R News 6 (2006). 

[51] S. Xiang, J. Compu t. Appl. M ath. 20 6, 688 (2007), ISSN 
0377-0427, URL |http : //portal . acm . org/citatioii7| 
cf m?id=1265610 . 1265788 

[52] G. A. Evans and J. R. Webster, J. Com put. Appl. Math 
112, 55 (1999), ISSN 0377-0427, URL |http : //portal . 

| acm . org/citation . cf m?id=335447 . 335464| 



